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ABSTRACT 

Experimental  and  theoretical  results  of  a nuclear  magnetic  double 
resonance  interferometric  study  of  a model  AX  spin  system  are  presented. 
Measurements  of  the  characteristic  relaxation  times  of  off-diagonal  density 
matrix  elements  corresonding  to  magnetic-dipole-forbidden  transitions  are 
presented,  and  the  use  of  such  relaxation  time  constants  to  obtain  information 
including  cross-correlations  of  the  fluctuating  fields  at  the  A and  X nuclear 
sites  is  discussed.  Off-resonance  effects  produce  large  changes  in  the  spectra, 
and  exploitation  of  these  effects  to  indirectly  measure  the  precise  resonance 
frequency  of  a spin  with  a small  magnetogyric  ratio  is  discussed. 


I.  INTRODUCTION 

' I 

This  paper  presents  a detailed  analysis  of  the  application  of  a 
recently  presented  interferometric  spectroscopic  technique^  to  a model  AX 
spin  system.  The  density  matrix  for  such  a coupled  spin  h system  is  a four- 

l 

by-four  matrix,  and  when  such  a spin  system  in  a strong  magnetic  field  is  j 

i 

prepared  in  a nonequilibrium  state,  the  subsequent  relaxation  of  the  system  to 
thermal  equilibrium  can  be  described  by  the  decay  of  the  elements  of  the  density 

i 

matrix,  where  the  diagonal  elements  relax  with  characteristic  time  constants, 

called  Tj's,  and  the  off-diagonal  elements  relax  with  characteristic  time  constants 

called  T2'sv  ’ . The  determination  of  all  of  the  Tj's  requires  a sequence  of 

experiments,  and  this  has  been  discussed  recently  in  detail  for  an  AX  system  by 

( 4 \ 

Mayne,  Adlerman,  and  Grant'  , while  measurement  of  the  off-diagonal  relaxation 

( 

\ 

rates  corresponding  to  magnetic-dipole-allowed  transitions  involves  application 

i 

of  a standard  Carr-Purcell  sequence.  This  leaves  off-diagonal  rates  corresponding 

. 

to  magnetic-dipole- forbidden  transitions  to  be  measured,  and  the  present  paper 
presents  a simple  method  for  direct  measurement  of  these  remaining  relaxation 
rates.  These  rates  are  of  particular  interest  since  they  can  contain  cross- 
correlation information  about  the  fluctuating  fields  at  the  sites  of  the  coupled 
A and  X nuclei.  In  our  review  of  the  literature  we  have  found  no  reports  of 
previous  efforts  to  measure  these  relaxation  rates  in  coupled  spin  h systems. 


although  we  do  want  to  call  attention  to  an  effort  to  detect  such  dipolar- 

forbidden  transitions  by  a multiple  step  excitation  process,  and  that  effort 

(5) 


27 

involved  a spin  5/2  nucleus,  A1  in  A1203V 

In  addition,  the  analysis  of  off-resonance  phase  effects  measured  in  the 
spectra  demonstrate  that  one  can  use  such  an  interferometric  scheme  to  measure 
indirectly  the  resonance  frequency  of  one  of  the  coupled  nuclei  by  observation  of 


ui 


3 


the  spectra  of  only  the  second  nucleus.  The  limiting  precision  of  such  an 
indirect  measurement  will  be  shown  to  be  equivalent  to  that  of  a direct  mea- 
surement. This  could  be  of  use  when  the  magnetogyric  ratios  of  the  two 
coupled  spins  differ  by  a large  amount, and  the  direct  observation  of  the  nucleus 
with  the  smaller  magnetogyric  ratio  is  made  difficult  by  poor  signal-to-noise 
conditions. 


4 


II.  DESCRIPTION  AND  EXPLANATION  OF  THE  EXPERIMENT 
A.  General  Description 

13  1 

The  system  we  chose  for  this  experiment  was  a conventional  AX(  C - H) 

13  1 

system  where  both  nuclei  are  spin  h-  The  C (S  spin)  and  H (I  spin)  were 

in  a liquid  and  were  weakly  coupled  to  each  other  via  a scalar  coupling  of  the 

form  J I • S.  Since  the  sample  was  in  a strong  magnetic  field,  oriented  along 

the  z-axis,  only  the  secular  part  of  the  scalar  interaction,  J Iz  Sz» contributed 

to  first  order  to  the  energy  level  spacings  for  this  system.  An  appropriate 

energy-level  diagram  is  shown  in  Figure  1.  The  four  levels  are  the  eigenstates 

of  the  Zeeman  Hamiltonian  corresponding  to  eta,  aB.  Ba,  and  BB,  where  the  a and 

B represent  the  spin  "up"  and  "down"  states  with  the  z-component  of  angular 

momentum  equal  to  hfi  and  - % ti,  and  the  first  Greek  letter  refers  to  the  I 

spin  while  the  second  refers  to  the  S spin.  Note  that  the  secular  part  of 

the  scalar  interaction  alters  the  energy  levels  slightly  (exaggerated 

considerably  in  Figure  1 for  purposes  of  explanation),  thus  giving  rise  to 

four  inequivalent  magnetic-dipol e-all  owed  transitions.  The  *H  NMR  spectrum 

is  a doublet  corresonding  to  transitions  1-3  and  2-4  (single  arrows),  and  the 
13 

C spectrum  is  also  a doublet  corresponding  to  transitions  l-2and  3-4  (double 

arrows).  This  doublet  splitting  is  crucial  for  the  present  experiment  since 

it  enables  us  to  selectively  irradiate  certain  transitions. 

The  experiment  is  shown  schematically  in  Figure  2.  First,  we  apply  a 

nonselective  tt/'2  pulse  (2  psec)  to  both  the  proton  transitions.  Then  immediately 

after  the  tt/2  pulse  we  apply  a low-power,  selective  pulse  of  length  x to  only 
13 

one  of  the  C transitions.  At  a time  aT  after  the  original  proton  tt/2  pulse 
we  apply  another  nonselective  tt  pulse  (4  psec)  to  both  proton  transitions,  causing 
a proton  spin  echo  to  be  formed  at  time  2AT.  The  decay  taken  from  time  t = 0 


(at  the  middle  of  the  echo)  is  recorded  for  Fourier  transformation  (FT).  The 

five  spectra  shown  in  Figure  3 are  such  FT  spectra  for  values  of  x = 0,  28,  56, 

84,  and  112  msec.  The  top  spectrum, for  which  t = 0, corresponds  to  a conventional 

spin  echo  pulsed  NMR-FT  proton  spectrum.  The  sample  was  a solution  of  9\% 

13 

C-enriched  sodium  formate  (NaCHO^)  dissolved  in  D^O  with  a small  amount  of 
1 113 

H impurity.  It  is  the  H and  the  C in  the  enriched  formate  ions  that  comprise 
our  AX  spin  system  and  give  rise  to  the  symmetric  proton  doublet  split  by 
J = 195  Hz.  The  small  peak  at  the  center  of  the  doublet  is  due  to  the  protons 

IP 

in  the  9 % of  the  formate  ions  containing  spinless  C.  The  large  peak  on  the 

far  left  is  due  to  the  small  amount  of  proton  impurity  in  the  D20.  When 

preparing  the  solution,  no  efforts  were  made  to  remove  dissolved  oxygen. 

When  looking  at  the  spectra  of  Figure  3,  two  effects  stand  out.  First, 

we  see  that  the  amplitude  of  the  formate  doublet  actually  goes  to  zero,  then 

negative,  back  to  zero,  and  finally  positive  again.  Secondly,  the  amplitudes 

of  both  members  of  the  doublet  do  not  remain  the  same.  The  former  effect  can 

be  explained  completely  by  the  spin  dynamics  of  the  system^  while  the  latter 

effect  can  be  explained  only  by  invoking  relaxation  phenomena.  Therefore, 

in  order  to  understand  the  results  of  this  experiment,  we  must  theoretically 

understand  the  spin  dynamics  of  such  a system  first  without,  and  later  with, 

relaxation  effects  taken  into  account. 

B.  Spin  Dynamics  Ignoring  Off-Resonance  and  Relaxation  Effects 

We  first  examine  the  spin  dynamics  of  our  experiment,  ignoring  relaxation 

effects  as  well  as  effects  of  the  *3C  selective  pulse  being  "off-resonance." 

There  are,  in  fact,  two  such  off-resonance  effects  to  be  ignored.  The  first 

is  the  possibility  of  the  selective  r.f.  radiation  being  at  a slightly  different 
13 

frequency  from  the  C transition  frequency.  The  second  is  the  possibility 


6 


that  although  the  selective  radiation  is  weak,  it  still  is  strong  enough  to 

11 

somewhat  disturb  the  other  C transition. 

Before  we  apply  the  first  pulse,  we  assume  that  our  spin  system  is  in 
thermal  equilibrium  and  that  it  can  be  described  by  the  following  4x4  density 
matrix: 


equil 


A 0 0 0 
0 B 0 0 
0 0 C 0 
0 0 0 D 


where  the  basis  states  are  chosen  to  be  the  states  1,  2,  3,  and  4 of  Figure  1 
(aa,  a$,  0a,  and  03,  respectively).  Note  that  the  off-diagonal  matrix  elements 
ai^  zero,  meaning  there  is  no  statistical  phase  coherence  in  the  system  and, 
furthermore,  in  the  high  temperature  approximation 


A - B = C - D 


1 13 

As  usual,  the  observable  H and  C magnetizations  are  proportional  to  the 
expectation  values  of  the  dimensionless  spin  operators,  I and  S.  The 
expectation  values  for  I and  S corresponding  to  the  various  observable 
magnetizations  are  related  to  the  density  matrix  elements  as  follows: 


= Re  Pi 


"Via = -Im  p13 


' <Ix>24  = p24 


’ <Iy>24  = "Im  p24 


<Iz>13  = ^pll  ' p33^  * <Iz>24  = ^p22  " p44* 


<s  >,o  = Re  P, 


• <Sx>34  = Re  p34 


<Sy>12  = -Imp12 


* <Sy>34  = "Im  p34 


<Sz>12  = ^pll  " p22^  ’ <Sz>34  = **^p33  - p44' 


'11  H22; 


Because  of  the  hermiticity  of  the  density  matrix.  Rep..  = Rep,,  and 

I J vJ  • 

Imp..  = -Imp...  The  subscripts  refer  to  particular  levels  between  which 

I J J I 

transitions  give  rise  to  observable  magnetization.  Because  of  the  scalar  split- 
ting, we  can  distinguish  between  the  transverse  magnetization  <1^3  and 
< Ix>24 » etc*  Each  corresponds  to  one  of  the  peaks  of  the  doublet.  However, 
since  we  do  not  know  the  sign  of  J,  we  do  not  know  which  magnetization 
corresponds  to  which  peak.  Fortunately,  this  is  no  serious  drawback,  and  we 
will  say  more  about  this  point  later  in  the  paper. 

Initially,  since  the  off-diagonal  elements  are  zero,  there  is  no 
transverse  magnetization.  The  purpose  of  the  proton  tt/2  is  then  to  create 
some  transverse  magnetization  and  phase  coherence.  After  the  ir/2  pulse, 

(assumed  to  be  along  the  x-axis  of  the  proton  rotating  frame),  we  have  the 
following  matrix,  making  use  of  Equation  (2)  and  ignoring  effects  due  to  the 
finite  width  of  the  pulse: 


p(  0)  = 


Js(A+C) 


^(B+D) 


-»si  (A-C) 


hi  (A-C) 


%(A+C) 


-^i (A-C)  0 


hi (A-C) 


MB+D) 


(4) 


8 


i 


i 
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The  Hamiltonian  acting  on  the  system  during  the  weak  C pulse  is 
(for  H = 1): 

H = Wj  Iz  + Sz  + J Iz  Sz  + Sx  cos  wt  (5) 

13 

where  we  have  chosen  the  pulse  to  be  along  the  x-axis  of  the  C rotating  frame. 

1 13 

The  frequencies  ojj  and  are  the  Larmor  frequencies  of  the  H and  C,  oj  is 

13 

the  frequency  of  the  C pulse,  is  its  amplitude,  and  J is  the  scalar 
coupling  in  radians  per  second.  We  now  analyze  this  by  transforming  to  a doubly 
rotating  frame  in  which  the  Hamiltonian  is  static,  by  choosing  to  go  to  the 
interaction  frame  of: 

Ho  = + wp  lz  (6) 

where  Wp  is  the  reference  frequency  of  the  phase  detector  we  used  to  detect 
the  proton  traverse  magnetization.  After  transforming  to  this  frame  we  have 
the  following  Hamiltonian 


H1  = Ao)Iz  + h J Sz  + J Iz  + ujj  Sx 


Aw  = ojj  - o)p  (7) 

a>5  - w = h J 


ignoring  the  counter-rotating  r.f.  component  as  usual.  The  density  matrix 
p is  also  taken  to  be  in  this  same  rotating  frame,  and  for  the  sake  of 
simplicity  we  will  not  use  any  special  notation  to  indicate  this.  Equation  (4) 
is  still  correct  because  we  have  assumed  the  lab  frame  and  the  interaction 
frame  to  be  coincident  before  the  -n/2  pulse,  and  we  are  justifiably  assuming 
the  pulse  is  short  enough  to  ignore  effects  during  the  pulse.  The  density 
matrix  evolves  in  time  according  to  the  liouville  equation: 


1 


The  first  kind  of  previously  mentioned  off-resonance  effect  was  ignored  when 
we  chose  w - = h J»  from  Equations  (7),  because  the  C radiation  is  being 

applied  right  "on"  one  of  the  resonances  of  the  doublet.  The  second  kind  of 
off-resonance  effect  is  ignored  by  assuming  that 

J » or  |H'n  - H 1 22  I >:>  H'l2  (10) 

This  approximation  means  that  the  r.f.  perturbation  is  so  weak  that  it  does 

not  mix  levels  1 and  2,  thus  we  set  tf'12  = = However,  levels  3 and  4 

13 

are  still  very  strongly  mixed  by  the  C pulse. 

Because  this  experiment  involves  observing  only  the  I transverse  magne- 
tization, Equations  (3)  show  that  it  is  sufficient  to  calculate  only  p^ 
and  p24  (or  p3l  and  p 42),  and  thus  we  shall  only  worry  about  determining  these 
matrix  elements.  Using  Equations  (8)  and  (9)  we  can  write  the  following 
differential  equations,  and  their  initial  conditions  using  Equation  (4): 


= -i  (HJ  + Aw)  P13  + i p14 
= -i  (y  + Aw)  P14  + i h Wj  P13 

= i (J5J  - Aw)  p24  + i % w1  p23 

= i (^J  - Aw)  p23  + i h Wj  p24 


Pj3  (0)  = P24  (0)  = H i (A  - C) 

Pi/i  (o)  ~ (o)  - ® 


Solving  these  equations,  we  can  determine  the  relevant  matrix  elements  at 

1 Q 

the  end  of  the  C pulse  of  length  t 


(t)  = H i (A  - C)  e 


-i  (JjJ  + Aw)x 


(t)  = Js  i (A  - C)  e 


i (JjJ  - Aw)t 


COS(%  Wjl) 

cos[h  Wji) 


The  time  evolution  of  the  density  matrix  from  the  time  t through  the 
it  proton  pulse  at  AT,  which  we  take  to  be  along  the  x-axis  of  the  proton 
rotating  frame,  and  finally  up  to  the  middle  of  the  echo  at  2AT , is  quite 
straightforward  to  calculate,  and  we  merely  state  the  pertinent  results: 


P3I  (2AT)  = e 


i(ijJ  + Aw)t 


P10  (t) 


-i(*sJ  - Aw)t 

p42  (2AT)  = e P24  U) 


t . 
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we  see  that  the  effect  of  our  selective  C pulse  is  to  cause  the  proton 
magnetization  to  oscillate  sinusoidally  with  a frequency  of  So  by  choosing 

the  length  of  x to  correspond  to  a normal  2tt  pulse  (co^x  = 2tt)  we  see  that 
effect  is  to  multiply  the  amplitude  of  the  proton  doublet  by  cos  (tt)  = -1. 

This  effect  is  the  manifestation  of  the  spinor  character  of  the  pseudo-two- 
level  system  composed  of  levels  3 and  4,  and  this  phenomenon  has  been  discussed 
in  detail  elsewhere^1 K The  spectra  in  Figure  3 were  taken  for  values  of  x 
corresponding  to  values  of  tOjX  = 0,  it,  2tt,  3tt,  and  4tt.  Note  that  after  a full 
4tt  rotation  the  phase  of  the  proton  magnetization  has  come  back  to  itself 
again,  which  is  further  consistent  with  spinor  behavior.  Thus,  we  can  see 
that  this  idealized  spin  dynamical  approach  can  at  least  explain  the 
oscillatory  behavior  in  Figure  3.  However,  we  need  to  understand  quantitatively 
the  amplitudes  of  both  members  of  the  doublet,  and  this  can  only  be  done  by 
appealing  to  relaxation  effects. 


By  examining  Equations (15) , one  notes  that  the  magnetization  linearly 
oscillates  rather  than  precesses.  By  looking  at  Equations  (11)  we  see  that 
the  matrix  elements  and  are  coupled,  as  well  as  P24  and  P23.  Solving 
for  p^  would  show  that  p14«sin  (^t);  therefore,  we  see  that  p13  and  p^ 
are  mutually  oscillating  and  when  one  has  a maximum,  the  other  is  zero, 
and  vice  versa.  According  to  the  definition  of  the  density  matrix: 


ptJ  ■ ci  cj  <16> 

where  C.  and  C-  are  the  complex  quantum  mechanical  amplitudes  for  the  system 

' J 

to  be  in  the  states  i and  j,  and  the  bar  represents  the  mixed  state  which 
results  from  taking  a statistical  ensemble  average  over  the  system.  Thus, 
we  see  that  when  p^  ^ 0,  statistically  our  system  is  in  a linear  combination 
of  states  1 and  4.  This,  in  turn,  means  that  there  is  phase  coherence  or 
"magnetization"  corresponding  to  the  forbidden  transition  1-4.  The  reason 
we  do  not  see  this  magnetization  is  twofold.  First  of  all,  since  it  is  magnetic- 
dipole-forbidden,  the  process  is  second  order  and  the  probability  of  the 
transition  is  correspondingly  low.  Secondly,  the  1-4  transition  frequency  is 
at  about  o)j  + and  we  would  have  to  make  our  phase  detection  reference 
frequency  Wp  closer  to  this  "double  quantum"  frequency  in  order  to  observe  this 
magnetization.  Similarly,  p£^  and  p^  form  a complementary  pair  of  oscillating 
variables,  and  033^0  implies  phase  coherence  between  the  levels  2 and  3,  which 
corresponds  to  a magnetic-dipole-forbidden  "flip-flop"  transition.  In 
anticipation  of  our  later  discussion  of  relaxation  effects,  we  mention  that  if 
the  off-diagonal  element  p^  relaxes  with  a different  rate  than  p^*  we 
would  expect  the  amplitude  of  the  1-3  peak  of  the  doublet  to  be  different  from 
the  prediction  of  Equation  (15).  This  is  because  for  t t 0,  the  system 
has  a probability  of  being  in  a linear  combination  of  states  1 and  4 as  well 


13 


as  1 and  3.  Thus,  the  ratio  of  <Iy(2AT)>  to  < Iy(2AT)>x  _0  could  be  either 
smaller  or  larger  than  that  predicted  by  Equation  (15),  depending  on  whether 
or  not  the  relaxation  rate  of  p14  is  faster  or  slower  than  the  relaxation  rate 
of  pjg-  Again,  similar  statements  can  be  made  about  the  rates  of  p 24  and  p2g. 
Thus,  by  attempting  to  quantitatively  fit  the  amplitudes  of  the  peaks  of  the 
formate  doublet  for  different  values  oft,  we  can  in  fact  determine  the  relaxation 
rates  of  the  density  matrix  elements  pl4  and  p^^.  In  the  next  section  we 
pursue  the  effects  of  being  "off-resonance"  on  the  C frequency,  while  still 
ignoring  relaxation.  We  will  see  that  these  "off-resonance"  effects  cannot 
explain  the  fact  that  the  doublet  does  not  remain  symmetric. 

C.  Spin  Dynamics  Including  Off-Resonance  Effects  but  Ignoring  Relaxation 
Effects 

As  mentioned  earlier,  there  are  two  important  "off-resonance"  effects 
13 

of  the  selective  C pulse  to  be  considered.  These  both  fall  naturally  out 

of  the  mathematics,  if  we  choose  the  proper  interaction  frame.  We  start  with 

the  same  Hamiltonians  as  those  in  Equations  (5)  and  (6);  however,  we  will 

13 

choose  our  frequency  of  C irradiation,  oj,  slightly  differently.  In  this 
case  Equations  (7)  must  be  modified  so  that  now  after  transformation  to  the 
interaction  frame,  we  have  the  following  remaining  Hamiltonian: 

H1  = Au  I2  + (y  + Aw)  + J Iz  Sz  + Sx 

Aw  = - a)  (17) 

u>5  - w = y + 6w 

13 

Therefore  <$w  is  just  the  difference  in  frequency  between  the  applied  C 
radiation  frequency  and  the  frequency  of  the  one  member  of  the  doublet  we  are 


14 


intending  to  irradiate,  and  thus  <5oj  measures  an  "off-resonance"  effect. 


Using  Equations  (17)  we  can  now  write  the  matrix  representing  H‘ 


Js(J+6w+Aa)) 


0 Js(Su)-Auj) 


^2°^  to+Aw) 


13r  H 


We  can  examine  the  other  type  of  C "off-resonance"  effect  by  not 


making  the  assumption  of  Equation  (10).  Therefore,  we  are  saying  our 


radiation  is  substantial  enough  to  somewhat  disturb  the  other  member  of  the 


3C  doublet.  Thus,  it  is  no  longer  necessary  to  assume  that  J » and 


we  are  not  so  far  "off-resonance"  from  the  other  1JC  transition  that  it  is 


irrelevant.  So  proceeding  with  this  in  mind,  we  can  use  Equation  (18)  and 


the  Liouville  Equation  (8)  to  get  the  following  set  of  coupled  differential 


equations: 


f>13  = -i  (^  + Auj)  Pj3  + i (p14  - P23) 


Pj4  * -i  (y  + Aw  + 6co)p14  + i >sw j (p13  - p24) 


P24  = i (*s  J ~ Aw)p24  + 1 (p23  ” Pj4) 


P23  = i (4  J - Aw  + 6o))p23  + i (p24  - p13) 


The  initial  conditions  for  p^,  p^4,  p^,  and  p„  are  identical  to  those  in 


Equations  (11)  . The  solution  of  this  set  of  equations  is  rather  tedious,  but 


the  important  matrix  elements  can  be  shov  to  be: 
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-iAoix 


Pl3(-0  = H(A-C)  e (i(a  cos  £i+x  + c cos  ft  x)  + (b  sin  ft+x  + d sin  ft  x)} 

-l/wr  (20> 

P24(T)  = *s(A-C)  e { i (a  cos  ft+x  + c cos  ft_x)  - (b  sin  ft+x  + d sin  ft  x)} 


where  we  have 


!±  = { GjO  + + <5oj) 2 + Wj2  ± t(^)  + ( j + <Sco)  + 00/ 


i 2 , 2 , 2 2 2 , 2 h 

- 2(j)  (75-  + 6w)  +2  (75-)  Wj  + 2 oj1  (j  + 6oj)  ] } 


2 % % 


where  ft+  refers  to  the  top  sign  and  ft_  refers  to  the  bottom  sign.  Also,  we 
have  the  following  values  for  a,  b,  c,  and  d in  Equations  (20): 

".2  - (f)? 


n_  - ft+ 

, 2 2 


i *■  c- 

(?)  - ft+ 
ft  2 - ft  2 


(f)m_2  - (f)2  - a)/] 

ft+  ( ft_2  - ft+2) 

2 2 2 

(?)  [(f)  + «*>!  - ] 

ft_  (ft_2  - ft+2) 


The  following  relations  also  hold  true 


a + c = l 


ft+b  + ft.d  = f 


16 


13 

After  the  C pulse  ends,  the  density  matrix  evolution  is  identical  to  that 
in  Equations (13).  Then  using  Equations  (13),  (20),  and  (3)  we  compute  the 
proton  transverse  magnetization  at  2AT  to  be 

<Iy(2AT)>13  = <Iy(2AT)>24  = <Iy(2AT)>T=0  {(b  sin  n+x  + d sin  fl_x) 

+ (a  cos  fi+x  + c cos  ft_x)  cos  (^  x)  } 

(24) 

< I x ( 2AT ) > j 3 = -<Ix(2AT)>24  = <Iy (2AT)>t=0  {(b  sin  n+x  + d sin  fi_x) 

- (a  cos  ft+x  + c cos  ft_x)  sin  (|-  x)  } 

where  <Iy (2AT)>t_q  is  again  the  amplitude  of  either  of  the  peaks  of  the 
doublet  for  the  case  where  x = 0.  Notice  that  the  y-components  of  the 
magnetization  of  the  proton  doublet  are  identical,  while  the  x-components 
merely  differ  by  a sign. 

Careful  examination  of  Equations (24)  tells  us  several  important  things. 

13 

First  of  all,  we  see  that  neither  of  the  C "off-resonance"  effects  we  have 
considered  can  possibly  lead  to  the  discrepancy  in  the  amplitudes  of  the  peaks 
in  the  doublet  of  the  spectra  in  Figure  3.  Any  such  off-resonance  effects 
may  alter  the  overall  amplitudes  of  the  peaks  as  well  as  introduce  dispersion 
to  the  peaks,  via  the  x-components,  but  the  two  peaks  must  remain  mirror 
images  of  each  other  in  the  spectrum.  Because  of  this,  we  can  then  rule  out 
magnetic  field  inhomogeneity  as  a possible  cause  of  that  discrepancy. 

Magnetic  field  inhomogeneity  can  be  treated  by  summing  a distribution  of 
peaks  being  off-resonance  by  different  amounts.  However,  any  distribution  of 
frequency  would  still  lead  to  overall  lineshapes  for  the  two  peaks  which  are 
mirror  images  of  each  other,  thus  ruling  this  effect  out  as  a possible 


explanation.  Also  we  can  precisely  determine  6w»  and  consequently,  the 

13 

resonance  frequency  of  the  C line  by  observing  the  effects  of  dispersion 
on  the  lineshape. 

The  full  behavior  of  the  functions  in  Equation  (24)  is  quite  involved. 

However,  we  can  look  at  certain  limiting  cases  to  separate  the  two  kinds 

of  off-resonance  effects.  To  try  to  understand  the  first  type  of  off-resonance 
13 

effect  where  the  C radiation  is  slight  "off"  the  intended  frequency,  we 

can  take  the  limit  of  Equations(24)  in  the  case  where  we  assume  the  approximation 

in  Equation  (10)  to  be  valid.  In  this  case  Equations  (24)  become: 

<Iy(2AT)>i3  = <Iy(2AT)>24  = <Iy(2AT)>T=Q  {cos  (h  jC\l  + 6a,2  t)  cos(^) 


+ 


+ 6a,2  x)  sin  (^)> 
(25) 


<Ix(2AT)>i3 


-<IX(2AT)>24 


2 x)  sin  (Sf-) 
) cos  (^p-) } 


For  the  case  where  6a,  « u^,  we  see  that  Equations  (25)  are  the  same  as  those 

of  Equations  (15)  except  that  there  is  a phase  error  of  6o>t/2  introduced  into 

the  peaks.  Thus,  both  peaks  have  the  same  amount  of  dispersion  mixed  in,  and 

the  sign  of  that  dispersion  is  opposite  for  the  two  peaks.  We  note  that  we 

13 

can  use  this  effect  to  our  advantage  in  determining  the  position  of  the  C 

13 

resonances.  We  could  arbitrarily  increase  t until  any  C off-resonance,  6o,» 
no  matter  how  small,  would  lead  to  noticeable  phase  v.nanges  in  the  proton  spectrum; 
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however,  an  upper  limit  on  x,  equal  to  approximately  the  inverse  of  the  13C 

1 3 

line  width,  effectively  limits  the  resolution  of  6co  to  the  C linewidth. 

Thus,  this  indirect  method  of  determining  the  position  of  the  C resonance 

seems  to  be  superior,  in  this  case,  to  other  methods  such  as  spin  tickling^6)  , 

13 

whose  resolution  is  ultimately  limited  by  the  proton  rather  than  the  C 
1 inewidth. 

The  second  type  of  off-resonance  effect,  in  which  we  take  into  account 

13 

the  effect  of  weakly  irradiating  the  other  member  of  the  C doublet,  can  be 
examined  by  assuming  that 

6uj  = 0 (26) 

but  not  making  the  assumption  of  Equations(lO) . Taking  the  subsequent 
limiting  case  of  Equations(24)  expanded  to  second  order  in  the  parameter 
x = ojj/J,  we  have  the  following 

<Iy(2AT)>13  = <Iy(2AT)>24  = <Iy(2AT)>T=0  {cos  (^x)  [cos [ (^x) (***) ] 
-hx2  sin  [Gjrr)(l  + %x2)]  sin  (^x)  ] 

+ x sin  (Hwjx)  sin  [(^x)(l  + *sx2)]  cos  (^x)  } 

(27) 

<Ix(2aT)>13  = -<IX(2AT)>24  = <Iy(2AT)>x=0  (cos  (s^x)  [sin[ (^jx) (»sx)] 

-*5X2  sin  [(^x)(l  + hx2)]  cos  (|x)  ] 

- x sin  (Jaojx)  sin  [ (^x) (1  + 4x2)]  sin  Gjx)  } 

Thus,  we  see  for  the  case  of  x « 1 (or  ^ « J)  that  Equations  (27)  are  the 
same  as  Equation  (15)  except  that  there  is  a phase  error  of ((dji/2)*  (x/2) 
introduced  into  the  peaks.  Thus,  again  both  peaks  have  equal  and  opposite 
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amounts  of  dispersion  mixed  in.  If  we  intend  to  use  the  first  type  of  off- 

13 

resonance  effect  to  find  the  C resonance  position  by  observing  the  amount 

of  dispersion  mixed  in,  we  must  take  into  account  also  this  other  phase  error 

due  to  the  second  type  of  off-resonance  effect.  This  does  not  hurt  our  resolution 
13 

of  the  C resonant  frequency,  but  it  does  mean  that  we  must  subtract  off  this 

second  effect.  Although  we  do  not  show  the  results  here,  computer  calculations 

of  Equations  (24)  indicate  that  if  both  types  of  off-resonance  effects  are 

present  simultaneously  and  they  are  both  small,  their  phase  errors  will,  in 

fact,  merely  add  linearly  rather  than  combine  in  some  more  complicated  manner. 

Thus,  it  seems  that  the  detailed  spin  dynamics  of  the  experiment,  taking 

off-resonance  effects  into  account,  seem  to  provide  the  necessary  insight  to 

13 

use  this  technique  as  a very  accurate,  indirect  method  of  determining  the  C 
resonant  frequency.  However,  only  in  the  next  section,  where  we  investigate 
the  effects  of  relaxation,  can  we  quantitatively  fit  the  measured  peak 
amplitude  to  theoretical  calculations. 

D.  Spin  Dynamics  Including  Relaxation  Effects  but  Ignoring  Off-Resonance 
Effects 

In  order  to  make  the  relaxation  calculation  more  tractable  and  more 

13 

easily  interpretable,  we  have  chosen  to  ignore  both  types  of  C off- 

121 

resonance  effects  in  this  section.  According  to  Redfield's  theory  of  relaxation'  ' , 
we  can  include  relaxation  effects  by  assuming  that  we  have  a modified 


Liouville  equation: 

it  pmn  = 1 £p’^' ^mn  ' T 


(m  ^n) 


Since  the  relaxation  times  are  long  compared  to  1/J,  we  can  assume  that  the 
damping  term  in  Equation  (28)  for  off-digonal  matrix  elements  depends  only 


A 


I 


» 
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on  the  same  matrix  element  appearing  on  the  left-hand  side  of  Equation  (28). 

Each  different  matrix  element,  pmn,  however,  is  assigned  its  own  rate  constant, 

1/T  . This  condition  of  Tmn  » 1/J  merely  means  that  we  do, in  fact,  have 

a well -resolved  proton  doublet  to  begin  with.  Furthermore,  the  requirement 

is  that  the  real  relaxation  time  Tm„  » 1/J  and  not  the  apparent  T * » 1/J 

mn  rr  mn 

★ 

where  Tmn  is  the  decay  constant  taking  magnetic  field  inhomogeneity  into 
account.  Thus,  this  criterion  is  easily  satisfied  in  our  experiment. 

Next,  by  using  Equations (28)  and  also  Equations  (11),  which  are  the  ones 
appropriate  for  ignoring  off-resonance  effects,  we  can  determine  the  fol lowing 
differential  equations: 


p13  = -1  (h  J + Aw)p13  + i ^ p14  - T p13 

Pj4  = -1  (^  J + Acl))pj4  + i p^3  - y^  Pj4 

f>24  = 1 (%  J - Ao))p24  + i P23  “ p24 

p23  = 1 (**  J " Au})p23  + i P24  ‘ T^3  p23 


(29) 


Again,  the  initial  conditions  are  the  same  as  in  Equations  (11).  The  relevant 
matrix  elements  can  then  be  determined  to  be 


-ify+AioH  + y^- )x 

P13(t)  = y(A-C)  e e 1 13  1 14 


{cos[Js  - y^-)2  T] 

1 13  14 


1_  - _1_ 

T14  T13 


/ 2 JL  / 

r (A)l  Vt f ) 

1 13  14 


(30) 
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p24^  = ^(A-C)  e 


IpsJ-Auh  M-t-  + — )T 

~ 1 0/1  * oo 


e ‘24  '23 


/ {T~'~T~  ] 

X {c.s[U,2-(Tl--Tl-)2TH  23  24 

i i24  i 23 


fi  tl  i \2 

^U1  'T  " T ) 

1 i 24  i 23 


sin  [h  / 2 - (— - - y^-)2  t] } 

i '24  '23 


As  in  the  other  cases,  the  Equations  (13)  for  the  evolution  of  the  density 
13 

matrix  after  the  C pulse  are  still  applicable.  Then  by  using  Equations 
(13),  (30),  and  (3),  we  get 


<Iy(2AT)>13  = <Iy(2AT)Vo  e-«T^  - {cost*  l*  - (tL  . J-)2  x 1 


<r--  r“> 

14  13 


14  '13 


/»? . (tl.  . ^ 

1 14  13 


sin[ 


, Tl  TT? 

4 /w,  - (= — ) 


2 ,U1  VT  " T 
1 14  13 


r in 

(31) 


<Iy(2AT)>24  = <Iy(2aT)>t=0  e'^-  T^)t  {cos ft  Z}  - - J-)2  t] 


23  '24 


'23 


r~> 

23  '24 


--  24  -sin  [h  A*  - (y^-  - =4-)2  t 
1 1x2  '23  ' 24 


1) 


<Ix^2AT)>13  = <Ix(2AT)>24  = 0 


where  we  note  that  <I^(2AT)>t_q  now  implicitly  includes  an  overall  relaxation 
term  exp  (-2AT/T13). 
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After  examination  of  Equations  (31),  we  see  that  the  effects  of  the 
relaxation  are  three-fold.  First,  we  see  that  there  is  an  exponential  damping 
or  increasing  of  the  amplitude  of  the  peak,  depending  on  the  difference  of 


the  two  relaxation  rates.  For  the  case  of  <IW(2AT)>-, , we  see  that  the  crucial 

13 


-y' 13 

quantity  is  (l/T..  - 1/T..J.  If  the  two  time  constants  T..  and  T,^  are 


2^  — x/  X 1 UIC  time  tuno  toil  to  '24 

equal,  then  we  see  no  effects  and  the  magnetization  is  the 
same  as  that  in  Equation  (15).  If,  however  < T13,  then  there  will  be  a 
damping  of  the  amplitude;  and  if  < T^,  then  there  will  be  an  increase 
in  the  peak's  amplitude.  The  second  effect  of  the  relaxation  is  to  alter 
the  oscillation  frequency;  however,  this  effect  is  second  order  in  the 
parameter  R13  = (1/Tj4  - 1/T13)/ojj.  The  third  effect  of  the  relaxation  is 
to  introduce  a first  order  (in  R13)  term  which  oscillates  like  sine  rather 
than  cosine.  The  main  effect  it  has  is  to  uniformly  translate  the  zero- 
crossings  of  the  amplitudes  as  it  effectively  introduces  an  overall  phase 
error  into  the  oscillating  term.  This  phase  of  the  oscillation  is  not  to  be 
confused  with  the  dispersion  - related  phase  discussed  extensively  in  Section  II  C 


of  this  paper.  Although  we  have  not  calculated  in  detail  the  effects  of 
13 


relaxation  and  x C off-resonance  effects  simultaneously,  we  feel  that  we 
can  safely  predict  the  results  for  small  off-resonance  parameters  by  a linear 
combination  of  the  dispersion  phase  errors  predicted  by  Equations  (25)  and 
(27)  and  the  magnetization  predicted  by  Equations  (31). 

Fitting  the  measured  amplitude  of  the  peak  to  the  theoretical  Equations  (31) 
can  yield  a value  for  R13,  and  knowing  idj  and  T13,  this  will 
then  yield  the  relaxation  time  T14.  A similar  procedure  for  the  other  peak 
will  subsequently  yield  T23  if  we  know  T^^.  It  seems  reasonable  to  assume  that 


T13  = T24 


(T2)i 


(32) 


23 


where  this  just  means  that  the  relaxation  time  constants  of  the  two  members  of 
the  proton  doublet  (corresponding  to  the  allowed  transitions  1-3  and  2-4)  are 
the  same.  Furthermore,  since  these  are  just  the  conventional  proton  T2's 
of  such  a system,  we  call  these  times  (T2)j. 

In  the  next  section  we  quantitatively  compare  theory  and  experiment 
and  determine  values  for  and  T23>  We  also  look  briefly  at  some  data 
which  show  off-resonance  effects  discussed  in  Section  II  C. 


III.  PRESENTATION  AND  DISCUSSION  OF  RESULTS 
A.  Relaxation  Data 

In  this  section  we  first  compare  the  spin  dynamical  theory  including 

relaxation,  of  Section  II  D,  to  experimental  data.  First,  we  need  to  rewrite 
Equations  (31)  in  a slightly  more  usable  form 


<Iy(2AT)>13  ^ Rn  

= e"'20R13  (cos (%6  /l-R  2)  + sin^e  /1-R  2)} 

<v24T)>t=o  13  13 


<Iy(ZAT)>24 

<Iv<24T)>I=0 


-Woil  f fl,r\  An  2\ 


e 2 24  {cos(^0  /l-R24  ) + 


sin(»se  /l-R24^)} 


where  we  have  now  normalized  the  amplitude  of  each  peak  to  the  value  for  the 
case  of  t = 0.  The  parameters  in  Equations  (33)  have  the  following  values 

0 = IdjT 


= (_1 = i_  /_1 

13  W1  T14  T13  “l  T14 

= i_  ( 1 L.)  = i_  . 

24  mj  T23  T24  (Oj  T23 


(f->  l) 

i2  1 


(f)j) 
'2  1 


24 


The  assumption  of  Equation  (32)  has  been  used  to  justify  the  substitutions 

for  T^3  and  ™ Equations  (34)  above. 

In  Figure  4 we  have  plotted  experimental  data  and  the  theoretical  Equations 

(33)  for  various  values  of  the  R parameters.  Figure  4A  shows  the  experimentally 

measured  amplitude  of  the  peak  of  the  proton  formate  doublet  which  was  not 

attenuated  for  long  values  of  t (the  right-most  member  of  the  doublet  in 

Figure  3).  Figure  4B  shows  the  experimental  amplitudes  of  the  other  peak  of 

the  doublet  which  was  strongly  attenuated  for  long  values  of  t (the  left-most 

member  of  the  doublet  in  Figure  3).  In  both  Figures  4A  and  B the  ratios  of 

the  amplitudes  of  the  doublet  peaks  to  the  amplitude  for  the  special  case  of 

t = 0 were  plotted  in  order  to  facilitate  direct  comparison  with  Equations  (33). 

Amplitudes  were  determined  from  areas  of  the  experimental  peaks,  and 

normalization  was  accomplished  by  comparison  of  the  1H  impurity  peak. 

The  triangles  indicate  data  points  for  experiments  where  the  echo  time,  2aT, 

was  chosen  to  be  90  msec  while  the  circles  indicate  data  points  for  experiments 

where  2aT  was  chosen  to  be  280  msec.  In  all  experiments,  the  same  value  for 

13 

the  strength  of  the  C r.f.  pulse,  uj,  was  used. 

The  question  now  arises  as  to  which  of  the  theoretical  Equations  (33) 
corresponds  to  which  of  the  doublet  peaks.  The  ambiguity  arises  because  we 
do  not  know  the  absolute  sign  of  J;  hence,  we  cannot  be  sure  which  of  the 
proton  transitions,  1-3  or  2-4,  is  the  lower  frequency  transition.  We  nave 
the  possibilities  of  irradiating  either  the  high  or  low  frequency  member 

1 O 

of  the  C doublet,  and  the  system  might  re. ct  with  attenuation  of  either 

the  high  or  low  frequency  member  of  the  proton  doublet.  Analysis  reveals  that 

we  can  unambiguously  assign  T^  and  T23,  but  we  cannot  determine  the  sign  of 

13 

0.  Let  us  consider  the  four  cases  of  irradiating  the  high  C (hi<j)  or 
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low  C (lo$)  transitions,  and  observing  attenuation  of  the  high  H (hij) 

or  low  *H  (lOj)  transitions.  It  turns  out  that  hi$  - hi j and  lo<.  - lOj  both 

predict  that  T ^ < (T2)j,  w^e  hij  ~ ^°I  ancl  ^°s  ’ h’l  Precl''ct  that 

T14  < (T2h-  *n  our  case>  we  experimentally  observed  only  hi,.  - hi  j and 

lo^  - 1 o j so  we  know  that  < (T2)j. 

By  the  discussion  of  the  preceding  paragraph,  we  now  can  determine 

values  for  R from  the  data  in  Figure  4 without  worrying  about  the  subscripts 

on  R.  Figure  4A  has  two  theoretical  curves  (Equation  (33)  ) drawn  for  R = 0 

(solid  line)  and  R = 0.02  (dashed  line).  From  these  curves  we  determine  that 

R = .01  + .01.  The  value  for  R and  the  error  limits  are  somewhat  subjectively 

determined  by  noting  that  the  curves  for  R = 0 and  for  R = .02  seem  to  nicely 

bound  the  scatter  in  the  data.  In  Figure  4B  we  have  drawn  theoretical  curves 

for  R = .20  (solid  line),  R = .18  (dashed  line),  and  R = .22  (dotted  line). 

The  curve  for  R = .20  fits  the  experimental  data  best,  and  the  curves  for 

R = . 18  and  R = .22  bound  the  scatter  in  the  data  rather  well.  So  from 

this  we  conclude  that  R = .20  ± .02  for  Figure  4B. 

We  determined  from  the  zero  crossing  of  the  experimental  data  in 

Figure  4A  to  be  w]/2tt  = 17.9  Hz.  This  procedure  is  convenient  since, for  this 

unattenuated  peak,  relaxation  effects  are  small  and  Equations  (33)  become 

equivalent  to  Equations( 15) . Another  option  to  determine  the  value  for 

would  be  to  measure  the  length  of  a 2n  pulse  by  actually  observing  the 
1 1 

C magnetization.  By  conventional  methods  we  measured  (T2)j  to  be  160  msec. 
Thus,  from  Figure  4A  the  value  of  R = .01  ± .01  means  that  T.^  = 140  ± 20  msec. 
The  value  of  R = .20  ± .02  from  Figure  4B  means  that  T ^ = 35  ± 3 msec. 

Althouqh  not  crucial  to  the  experiment,  for  general  information  we  measured 
the  (T^)  (this  is  the  13C  T2)  to  be  480  msec,  the  (Tj)j  (normal  proton  T^) 
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to  be  4.2  sec,  and  (Tj)$  (normal  C T^)  to  be  8.7  sec.  All  of  these  times 
are  shorter  than  would  be  expected  in  a system  purged  of  dissolved  oxygen. 

Using  this  method,  it  is  now  possible  to  directly  measure  all  of  the 
relaxation  rates  pertaining  to  the  off-diagonal  density  matrix  elements  of 
such  a system.  In  our  case  the  only  four  unique  non-zero  rates  have  the  time 
constants  (T2  ^ I * S’  T14’  ancl  T23‘  The  Afield  theory^  says  that 
the  general  relaxation  of  the  density  matrix  is  of  the  form 

dt  pmn  Rmnij  pij 

1 J 

where  the  Rmnij  is  a super-matrix.  Thus,  the  relaxation  of  a particular 

element  of  the  density  matrix  in  general  depends  on  the  values  of  other  elements 

as  well.  However,  according  to  our  assumption  of  Equation  (28),  the  relaxation 

of  each  off-diagonal  element  depends  only  on  its  own  value.  This  is  not  the 

case  for  diagonal  elements  even  if  we  still  assume  that  the  relaxation  times 

( 4 1 

are  long  compared  to  1/J.  A recent  paper  of  Mayne,  et  al.  ' shows  a very 
pretty  experiment,  performed  on  the  same  chemical  system  we  used,  which  measures 
all  of  the  elements  of  the  super-relaxation  rate  matrix,  which  pertain 

to  diagonal  elements  of  the  density  matrix.  So,  by  performing  experiments 
like  these  as  well  as  the  experiments  described  here,  it  is  possible  to 
map  out  the  values  of  the  entire  16  x 16  relaxation  matrix  for  this  AX  system. 
Generalized  versions  of  this  phase  interferometric  spin  spectroscopic  technique 
are  easy  to  generate  for  more  complicated  systems  such  as  AX^,  etc. 

If  the  decay  times  of  the  peaks  had  been  long  compared  to  the  length  of 

1 O 

our  C pulse  (as  might  have  been  the  case  had  we  removed  dissolved  oxygen), 
there  would  have  been  several  alternatives.  First  of  all,  we  could  have  just 
increased  the  time  t until  relaxation  effects  were  noticeable.  Care  should 

i 
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be  used  in  choosing  the  proper  cj.  because  if  during  the  time  x the  value 
of  cos  (gjjt/2)  undergoes  many  cycles,  then  the  inhomogeneity  of  the  r.f. 
magnetic  field  (coj/yj)  could  cause  the  proton  magnetization  to  dephase,  thus 
giving  rise  to  an  effective  decay  rate  which  could  interfere  with  the  measurement 
of  the  relaxation  decay.  This  will  not  be  a problem  if  the  following  criterion 
is  maintained 


Awi 

— « R 
U1 


(36) 


where  R is  the  dimensionless  parameter  of  Equations (34) . The  symbol  Awj 

stands  for  the  variation  of  the  C r.f.  perturbation  across  the  sample,  so 

Acjj/ojj  is  the  fractional  r.f.  inhomogeneity,  determined  by  the  geometry  of 

the  r.f.  coil.  A second  alternative  to  the  initial  problem  would  be  to 

stop  the  C pulse  at  a time  when  the  proton  magnetization  has  completely 

disappeared.  This  would  correspond  to  values  of  e = WjT  = it,  3tt,  5tt,  etc. 

At  that  time  the  state  of  the  system  would  have  and  as  the  only  non-zero 

off-diagonal  matrix  elements,  and  the  system  would  relax  with  the  corresponding 

times  T^  and  T^^-  After  waiting  some  appropriately  long  time  one  would 

"retrieve"  the  proton  magnetization  by  the  appl ication  of  another  selective 
13 

C pulse.  Measurement  of  the  magnetization  would  then  yield  the  desired 

relaxation  times.  A method  similar  to  this  last  suggestion  was  used  in  a 

(5) 

quadrupolar  system  by  Hatanaka,  et  al.  ' to  measure  off-diagonal  forbidden 

27 

transition  relaxation  rates  for  spin  5/2  nuclei,  Al  in  Al^O^.  Of  these  two 
methods  the  former  should  work  in  all  cases,  while  the  later  is  applicable 


only  when 


for  i , j = 1,4  and  2,3 


(37) 
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B.  Qualitative  Off-Resonance  Data 

In  this  section  we  very  briefly  present  two  spectra  showing  some  of  the 

off-resonance  effects  of  the  spin  dynamics  described  in  Section  II,  C. 

13 

Figure  5 contains  spectra  taken  for  ojjT  = 2tt  in  which  the  C radiation  is 

slightly  off-resonance.  Both  types  of  off-resonance  effects  contribute  to  these 

spectra.  For  the  value  of  stated  earlier,  we  have  the  important  parameter 

of  Equation  (27),  x = .092.  This  means  that  if  the  C offset  frequency, 

Sw,  is  zero,  there  will  be  an  apparent  phase  error  of  about  8°  due  to  the 

second  type  of  off-resonance  effect.  The  first  kind  of  off-resonance  effect 

13 

then  predicts  that  irradiating  the  C line  at  6(o/2tt  = .8  Hz  should  produce 

spectra  with  no  dispersion  present.  This  can  be  estimated  by  setting  the 

apparent  phase  errors  from  Equations  (25)  and  (27)  equal  to  each  other  to  get 
2 

oj, 

6u  = -jj-  (38) 

13 

Thus,  we  speak  of  this  C irradiation  frequency  as  being  the  apparent 
resonance  frequency. 

13 

Figures  5A  and  B show  spectra  taken  for  the  values  of  the  C 

irradiation  frequency  above  and  below,  respectively,  the  apparent  resonance 

frequency.  These  spectra  were  taken  10  Hz  apart,  so  we  would  expect  to  see 

a phase  difference  of  k <Wr  equal  to  about  101°.  It  is  difficult  to  determine 

from  the  spectra  exactly  what  the  phase  errors  are  because  this  was  not  done 

with  a high  resolution  spectrometer,  and  thus  the  dispersion  phase  signals  can 

overlap  somewhat.  However,  we  can  still  estimate  that  Figure  5A  is  about  4 Hz 

above  the  apparent  resonance  frequency  and  Figure  5B  is  about  6 Hz  below, 

corresponding  to  phase  errors  of  +40°  and  -60°.  Very  accurate  measurements 

13 

could,  in  principle,  be  made  to  very  precisely  determine  the  apparent  C 
resonance  frequency. 
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IV.  IMPLICATIONS  OF  THE  EXPERIMENT 

We  discuss  here  two  important  applications  of  the  experiment  described 

in  the  foregoing  sections.  One  of  these  is  the  precise  determination  of 

13  1 

C resonance  frequencies  through  the  observation  of  changes  in  the  H spectrum, 

as  described  in  Section  II,  C.  Since  the  sensitivity  of  the  technique  does 
not  depend  on  the  magnetogyric  ratio,  y<.,  of  the  S spins  it  potentially  is  a 
highly  desirable  method  for  the  detection  of  low-y  spin  resonances.  The  only 
condition  for  being  able  to  perform  the  experiment  is  that  the  doublet  of  the 
S-spectrum  is  well  resolved, and  under  the  reasonable  assumption  that  0 as  well 
as  the  linewidth  (due  to  field  inhomogeneity)  is  proportional  to  ys,  it  is 
clear  that  the  magnitude  of  y^  does  not  impose  any  limitation  on  the 
applicability  of  the  method.  The  other  main  application  is  the  measurement 
of  T^'s  associated  with  forbidden  transitions,  in  our  case  T23  and  T14. 

Although  we  will  not  give  an  extensive  treatment  of  the  various  relaxation 
mechanisms  that  possibly  might  prevail,  we  wish  to  illustrate  the  particular 
significance  of  these  relaxation  times  by  the  following  very  simple  model. 
Suppose  that  the  spin  relaxation  of  an  AX  system  is  caused  by  randomly  fluctu- 
ating local  fields  AHj  and  AHS  at  the  sites  of  the  nuclei  I and  S,  respectively, 
and  let  us  assume  that  these  fields  are  parallel  to  the  external  Zeeman  field. 
Defining  Awj(t)  = Yj  AH^,(t)  and  Aus(t)  = y$  AHs(t),  we  then  have  for  the 
fluctuating  random  Hamiltonian: 


r-’ 

mn 


where  the  bar  denotes  the  ensemble  average  and  x is  the  correlation  time 
of  the  fluctuations.  Hence,  we  find  for  this  simple  model 


/I  x _ 1 1 . 2 

VT“)  T " T T Aa)I  T 

i2  i i13  i24  i 


(j  )$  - T T “ T 

I2  b I12  1 34 


— = (Au)j^  + Aw^  - 2AcojAa)<.)  t 
23 


1 2 2 
j — = (Aa>j  + A + 2Aa)jAco^)  r 

While  the  normal  relaxation  times  depend  on  Auj2  and  Ac o^\  which  measure  the 
strengths  of  the  fluctuating  fields,  the  times  T23  and  are  in  addition 
related  to  Acoj  Acu^.  Thus,  measurement  of  the  T2's  of  forbidden  transitions 
provides  information  about  the  cross-correlation  between  the  two  local  fields. 
One  should  be  aware  that  this  result  strongly  depends  on  the  particular 
relaxation  mechanism  chosen.  For  instance,  fluctuating  fields  perperdicular 
to  the  Zeemen  field  with  very  short  correlation  times  give  rise  to  uniform 
relaxation  rates  for  all  the  off-diagonal  elements  of  the  density  matrix, 


= (tt)s  = t: 


23  '14 


= ^(Au).^  + Awc^)  T 

I 1 A 1 


(42) 
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Nevertheless,  we  feel  that  the  physical  picture  emerging  from  the  first 
example  will  have  general  implications,  and  that  knowledge  of  relaxation  times 
like  T23  and  will  be  of  help  in  the  determination  of  the  detailed  nature 
of  molecular  motions,  such  as  anisotropic  tumbling. 

Furthermore,  the  phase  interferometric  technique  is  in  general 
applicable  whenever  we  have  a system  with  two  or  more  inequivalent  transitions 
having  one  quantum  mechanical  level  in  common.  It  can,  for  .instance,  be 
applied  to  systems  consisting  of  two  coupled  nuclear  spins,  a nucleus  and  a 
free  electron, or  a nucleus  with  quadrupolar  interaction. 
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FIGURE  CAPTIONS 

Figure  1.  Energy  level  diagram  for  an  AX  spin  system.  The  a and  3 represent 

the  two  eigenstates  spin  up  and  spin  down  of  the  spin  h particle. 

The  first  Greek  letter  represents  the  state  of  the  *H  spin  and  the 

13  1 

second  represents  the  state  of  the  C spin,  so  that  two  H transitions 

13 

are  shown  with  single  arrows,  while  the  two  C transitions  are  shown 

with  double  arrows.  The  numbers  1,  2,  3,  4 are  used  to  refer  to 

the  various  energy  levels  or  to  the  eigenstates  to  which  they 

correspond.  The  relative  Zeeman  energies  for  *H  (56.4  MHz)  and 
13 

C (14.2  Mhz)  have  been  drawn  to  scale,  but  the  effects  of  the  weak 
coupling  have  been  greatly  exaggerated  for  emphasis. 

Figure  2.  Radio  frequency  pulse  sequence  used.  A tt/2  pulse  is  applied  to  both 

*H  transitions.  Then  a selective  pulse  of  length  x is  applied  to  only 

13  l 

one  of  the  C transitions.  A n pulse  is  then  applied  to  both  H 

transitions  at  a time  AT  after  the  tt/2  pulse,  which  creates  a spin 

echo  at  a time  2AT  after  the  tt/2  pulse.  The  signal  is  recored 

from  2AT(t  = 0),  defined  to  be  the  middle  of  the  echo  for  Fourier 

transformation. 

Figure  3.  Proton  phase  interferometric  spectra  for  different  values  of  x. 

The  doublet  split  by  195  Hz  is  due  to  the  *H  coupled  to  the  in 

13 

those  formate  ions  containing  C.  The  small  peak  at  the  center  of  the 

1 12 

doublet  is  due  to  the  H in  the  formate  ions  containing  spinless  C, 

while  the  large  peak  on  the  far  left  is  due  to  the  small  amount  of 

*H  impurity  in  the  solvent.  The  five  spectra  are  for  values  of 

x = 0,  28,  56,  84,  and  112  msec,  corresponding  to  values  of 

wjx  = 0°,  180°,  360°,  540°,  and  720°. 
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Figure  4.  Comparison  of  experimental  and  theoretical  values  of  the  amplitudes  of 
the  peaks  of  the  proton  doublet  vs.  the  value  of  e = co^t.  The 
amplitudes  of  the  peaks  are  plotted  on  the  vertical  axis,  normalized  to 
the  peak  amplitude  for  the  special  case  of  T = 0.  The  value  of 
( in  degrees)  is  plotted  on  the  horizontal  axis.  The  triangles  indicate 
experimental  data  for  which  the  echo  time  (2aT)  was  90  msec,  and  the 
circles  are  for  2aT  = 280  msec.  In  all  cases  the  same  value  of 
= 112  rad/sec  was  used. 

A.  This  solid  line  corresponds  to  the  theory  for  the  value  of  the 
parameter  R = 0.  The  dashed  line  represents  the  theory  for 

R = .02. 

B.  The  solid  line  represents  theory  for  R = .20,  the  dashed  line  is 
R = .18,  and  the  dotted  line  is  R = .22. 

Figure  5.  Proton  phase  interferometric  spectra  (for  0 = oj^t  = 360°)  showing 
13 

C off-resonance  effects. 

13 

A.  This  spectrum  was  taken  with  the  C irradiation  about  4 Hz 

13 

above  the  apparent  C resonance,  thus  the  doublet  shows  a phase 
error  of  about  + 40°. 

13 

B.  This  spectrum  was  taken  with  the  C irradiation  about  6 Hz  below 

13 

the  apparent  C resonance,  thus  the  doublet  shows  a phase  error 
of  about  - 60°. 
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